Black-Hole Bombs and Photon-Mass Bounds 
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Generic extensions of the standard model predict the existence of ultralight bosonic degrees of 
freedom. Several ongoing experiments are aimed at detecting these particles or constraining their 
mass range. Here we show that massive vector fields around rotating black holes can give rise to a 
strong superradiant instability which extracts angular momentum from the hole. The observation 
of supermassive spinning black holes imposes limits on this mechanism. We show that current 
supermassive black hole spin estimates provide the tightest upper limits on the mass of the photon 
(m v < 4 x 10~ 20 eV according to our most conservative estimate), and that spin measurements 
for the largest known supermassive black holes could further lower this bound to m v < 10~ 22 eV. 
Our analysis relies on a novel framework to study perturbations of rotating Kerr black holes in the 
slow-rotation regime, that we developed up to second order in rotation, and that can be extended 
to other spacetime metrics and other theories. 
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Introduction. The properties of matter making up our 
universe are mostly unknown. Strong evidence (e.g. from 
galactic rotation curves and from gravitational lensing) 
points to the existence of elusive, weakly-interacting mat- 
ter as the most abundant clement in the universe. An in- 
teresting possibility is the existence of ultralight bosonic 
degrees of freedom, such as those appearing in the "string 
axiverse" scenario P, 0, or of massive hidden U(l) vec- 
tor fields, that are also a generic feature of extensions of 
the standard model @-|(|. 

Massive fields around rotating black holes (BHs) can 
trigger a superradiant instability, the so-called "black 
hole bomb" [3]. This instability is well understood in 
the case of massive scalar fields [8Ml5| : it requires the ex- 
istence of negative energy states in a region around the 
BH known as the ergoregion. The instability is regulated 
by the dimcnsionlcss parameter Mu (from now on we set 
G = c = 1), where M is the BH mass and m s = uh 
is the scalar field mass, and it is most effective when 
Mu ~ 1 and for maximally spinning BHs. For a solar 
mass BH and a field of mass m s ~ 1 eV the parameter 
Mu ~ 10 10 , and therefore in many cases of astrophysical 
interest the instability timescale is larger than the age 
of the universe. Superradiant instabilities strong enough 
to be obscrvationally relevant (M u ~ 1) can occur cither 
for light primordial BHs that may have been produced in 
the early universe or for ultralight exotic particles 

found in some extensions of the standard model [H, . In 
the string axiverse scenario, massive scalar fields with 



10~ 33 eV < m s < 10~ 18 eV could play a key role in 
cosmological models. Superradiant instabilities may al- 
low us to probe the existence of ultralight bosonic fields 
by producing gaps in the mass-spin BH Rcgge spectrum 
[lj, l2[ , by modifying the inspiral dynamics of compact bi- 
naries [ill, [l^, H(| or by inducing a "bosenova" , i.e. a 
collapse of the axion cloud |2lH23j . 

The curved spacetime dynamics of massive vector fields 
has only been studied in nonrotating backgrounds (23 - 
27} • While superradiant instabilities are expected to oc- 
cur also for massive vector fields, quantitative investiga- 
tions have been hampered by our inability to fully un- 
derstand the (massive vector) Proca equation 



V CT F ff " - u 2 A p = , 



(1) 



where A M is the vector potential, F^ v = d^A v — d v A 
and m v 



the Lorenz condition V ^A^ 



ah is the mass of the vector field. Note that 
is automatically satis- 
fied and the Proca field propagates 3 degrees of free- 
dom (26[. In a nutshell, the problem is that Eq. ([1]) docs 
not seem to be separable in the Kerr background. 
Framework. The Proca perturbation problem in the 
Kerr metric becomes tractable if we work in the slow- 
rotation approximation. Let us focus on the Kerr metric 
in Boyer-Lindquist coordinates 
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and the angular momentum of the BH, respectively. In 
what follows, we shall expand the metric and all other 
quantities of interest to second order in 5. The proce- 
dure to separate the linearized perturbation equations to 
first order in a was first proposed by Kojima in the con- 
text of stellar perturbations |28l - l30l |. but it can be gen- 
eralized to any order in a and to generic (scalar, vector, 
tensor, etc.) perturbations of stationary and axisymmet- 
ric spacetimes. The details of the procedure will appear 
elsewhere [3l|; here we only present the main results. 

In the slow-rotation limit the perturbation equations, 
expanded in spherical harmonics and Fourier trans- 
formed in time, yield a coupled system of ODEs. In the 
case of a spherically symmetric background, perturba- 
tions with different harmonic indices (£, to), as well as 
perturbations with opposite parity, are decoupled. In 
a rotating, axially symmetric background, perturbations 
with different values of the azimuthal number m are still 
decoupled, but perturbations with different values of £ 
are not. However, in the limit of slow rotation there is a 
Laporte-like "selection rule" [32j ]: at first order in a, per- 
turbations with a given value of £ are only coupled with 
those with I ± 1 and opposite parity, similarly to the case 
of rotating stars. At second order, perturbations with a 
given value of £ are also coupled with those with £ ± 2 
and same parity, and so on. 

In general, the perturbation equations can always be 
written in the form [3l[ 

= At + amAe + a 2 At 
+ a{QiPi-i + Q e+ iV e+ i) 

+ a 2 [Qe-iQiAi-2 + Qt+2Qe+iAt +2 \ + 0(a 3 ) , (2) 

= Vi + amVi + a 2r Pt 
+ ~a{Q t Ai-x + Qe+iAe+i) 

+ a 2 \Qt-iQeVt-2 + Q^Qi+ifw] + 0{a 3 ) , (3) 



where Qg = y e - e2 ™* and the coefficients At and Vi (with 

various superscripts) are linear combinations of axial and 
polar perturbation variables, respectively. 

The general method can be specialized to the Proca 
equation (fTJ). We expand the vector potential as (2(| 



8A li (t,r,0,<p) = Y, 











+ 


.«(4)S|/A_ 





i)Y'/r 



'■(1) 



u\ 2) Y"/{rf) 
Y(/A 



'(3) 



where b = (<&, if), A = £(£ + 1), / = A/(r 2 + h 2 M 2 ), Y e 
are scalar spherical harmonics, = (Y#, Y~) and S e b = 
(Sq,SV) are vector spherical harmonics [3 if . whereas 
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(r, t) are polar 



and axial perturbations, respectively. Separating the an- 
gular variables, we find that Proca perturbations in the 
slow-rotation limit, up to second order, are described by 



two sets of equations [311 ] : 

£>A*A* + V A *A* = 0, (4) 

PpV+VpV = 0, (5) 

where £>a,p are second order differential operators, Va,p 
are matrices, M/a = ( w (4)> u (2) ' u (3) 



u^tf&Sugf) and 



(«| 2 j,Mpj,Uj 4 j t u (2) ' u (3) )• The function can be ob- 
tained from the Lorenz condition once the three dynam- 
ical degrees of freedom are known [3l|. When a = 0, 
the equations above reduce to Proca perturbations of 
a Schwarzschild BH (2(|. However, rotation introduces 
mixing between perturbations of different parity and dif- 
ferent multipolar indices. 

Numerical Results. Once suitable boundary condi- 
tions and a time dependence of the form e~ 1LJt are im- 
posed, Eqs. ((H)-© form an eigenvalue problem for the 
complex frequency uj = u)r + Physically moti- 

vated boundary conditions correspond to either quasi- 
normal modes (perturbations having ingoing-wave con- 
ditions at the horizon and outgoing-wave conditions at 
infinity [33j ) or bound states (perturbations that are spa- 
tially localized within the vicinity of the BH and de- 
cay exponentially at infinity). Here we focus on bound 
modes. By analogy with the scalar field case we would 
expect these modes to become superradiantly unstable 
for ujfj < mVtfj [To| . where il^ = <V(2r + ). Our analysis 
shows, for the first time, that massive vector fields do 
indeed become unstable in this regime. 

The bound state modes of the system (UJ-© can be 
found by standard numerical methods [3l| . When m > 
we find that, within numerical errors, the imaginary part 
of the modes has a zero crossing when 



(G) 



which corresponds to the onset of the supcrradiant 
regime. When ojr < jtiSIjj the modes are unstable and 
t = ljJ 1 is the instability growth timescale. Note that, 
although the supcrradiant condition (|6]) appears as a first 
order effect, in fact ujM ~ a at the onset of supcrradi- 
ance. The field equations contain terms proportional to 
auj and to uj 2 , so a second-order expansion is needed for 
a self-consistent study of the unstable regime [3l| . 

When Mu < 0.1 our data for the fundamental modes 
are consistent with a hydrogenic spectrum, ujr ^ [i and 



M loi ~ 7s£ (dm — 2r+ /i) (Mfi) 



U+5+2S 



(7) 



where 75^ is a coefficient that depends on £ and on the 
"polarization" index S, with S = for axial modes and 
S = ±1 for two classes of polar modes [26|. In the axial 
case the numerical results are also supported by an ana- 
lytical formula, that can be found by appl ying Starobin- 
ski's method of matching asvmptotic s IIOL |34| to Eq. (|4| 
at first order. A detailed calculation [31j yields precisely 
Eq. © with S = and 701 = 1/12. This is consistent 
with our numerical data, that yield 701 ~ 0.09 ± 0.03 
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(here and in the following, numerical errors are estimated 
by comparing the results at first and second order and by 
taking the maximum deviation between the fit and the 
data) . The instability timescale in the axial Proca case is 
four times shorter than in the scalar case. Similar results 
can be derived also for axial modes with £ > 1 and for 
the overtones. 

To our knowledge, this is the first estimate of the in- 
stability timescale of massive vector fields around spin- 
ning BHs. Although Eq. (J7J is strictly valid only when 
a <C 1 and M/i <C 1, in the case of massive scalar fields 
it provides estimates in good agreement with exact re- 
sults [H, [H up to a < 0.99: for example, Eq. © over- 
estimates the exact result by only 3% when a = 0.7, and 
by less than 70% when a = 0.99 [H|. 

In the massive vector case, the good agreement be- 
tween the first- and second-order calculations suggests 
that the slow-rotation expansion can be trusted even for 
moderately large spins, a < 0.7. Therefore it is reason- 
able to expect that extrapolations of Eq. ([7]) from the 
slow-rotation limit should at least provide the correct or- 
der of magnitude (and possibly a reliable quantitative 
estimate) of the instability timescale far from extremal- 
ity If we extrapolate Eq. (J7J) to a — > 1, we expect to 
overestimate the instability by about one order of mag- 
nitude [H|. 

As we shall discuss below, astrophysical bounds on 
vector field masses (scaling as -y~ 1 /( 4 ^+5+2S) [3lj] ) de- 
pend very mildly on uncertainties in the 757 coefficients. 
However they are sensitive to the scaling with \x in the 
M/j, <C 1 regime, and it is crucial to obtain reliable results 
in this limit. Unfortunately the calculation of unstable 
modes when M\i < 0.02 is challenging due to numerical 
inaccuracies, but the consistency between our numerical 
data in the axial case and the analytical formula (which 
is valid when M/j, <gc 1) is reassuring. The equations 
governing polar modes are much more complex, and we 
could not find analytical results supporting the fit 
when 5^0 (cf. [26[). Our data for S — — 1 polar modes 
at the onset of the i = m = 1 instability are consistent 
with Eq. © with 7_n ~ 20 ± 10, but some experimen- 
tation showed that different fitting functions can provide 
even better fits. Due to these uncertainties in the fit ([7]) 
for the polar case, in the following we shall discuss the 
consequences of our results mainly for axial modes, for 
which we derived the instability timescale also analyti- 
cally |3lj |. This choice is very conservative because, ac- 
cording to Eq. (JT)), polar modes with S = — 1 exhibit the 
strongest instability. Indeed, for fixed values of a and /1 
the instability of polar modes with S = — 1 is typically 
two or three orders of magnitude stronger than in the 
axial case [3lj . 

Astrophysical bounds on the photon mass. Our re- 
sults, together with reliable supermassive BH spin mea- 
surements, can be used to impose stringent constraints on 
the allowed mass range of massive vector fields. These 
bounds follow from the requirement that astrophysical 
spinning BHs should be stable, in the sense that the in- 



stability timescale r should be larger than some obser- 
vational threshold. For isolated BHs we can take the 
observational threshold to be the age of the Universe, 
THubbie = 1-38 x 10 10 yr. However, for supermassive BHs 
we may worry about possible spin growth due to merg- 
ers with other BHs and/or accretion. The most likely 
mechanism to produce fast-spinning BHs is prolonged 
accretion (35| . Therefore, a conservative assumption to 
estimate the astrophysical consequences of the instability 
is to compare the superradiance timescale to the (mini- 
mum) timescale over which accretion could spin up the 
BH. Thin-disk accretion can increase the BH spin from 
a = to a ~ 1 with a corresponding mass increase by 
a factor y/6 [36||. If we assume that mass growth oc- 
curs via accretion at the Eddington limit, so that the BH 
mass grows exponentially with e-folding time given by 
the Salpeter timescale rsaipctor = 4.5 x 10 yr, then the 
minimum timescale for the BH spin to grow via thin-disk 
accretion is comparable to rsaipctcr- 




M/M Q 



FIG. 1. Contour plots in the BH Regge plane 0| correspond- 
ing to an instability timescale shorter than rs a i pc tcr for dif- 
ferent values of the vector field mass m v — \iK and for axial 
modes with I — m = 1. Dashed lines bracket our estimated 
numerical errors, 701 « 0.09 ± 0.03 in Eq. ©. The experi- 
mental points (with error bars) refer to the supermassive BHs 
listed in Table 2 of [37| and the rightmost point corresponds 
to the supermassive BH in Fairall 9 [3SJ. Supermassive BHs 
lying above each of these curves would be unstable on an ob- 
servable timescale, and therefore each point rules out a range 
of Proca field masses. 

Brenneman et al. [13 have recently presented a list of 
eight supermassive BH spin estimates. In order to quan- 
tify the dependence of Proca field mass bounds on the 
mass and spin of supermassive BHs, in Fig.[TJwe show ex- 
clusion regions in the "BH Regge plane" (cf. Fig. 3 of 0). 
To be more specific, we plot contours corresponding to an 
instability timescale of the order of the Salpeter time for 
three different masses of the Proca field (m v = 10 -18 eV, 
2 x 10" 19 eV, 4 x 10" 20 eV) and for axial modes (S = 0). 
The plot shows that observations of supermassive BHs 
with 1O 6 M < M < 1O 9 M and a > 0.3 would exclude 
a wide range of vector field masses. 
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FIG. 2. Comparison between axial and polar I = m = 1 
instability windows for m v = 10 _20 eV. The right boundary 
of the instability window does not depend on uncertainties 
in the fits and it is given by J/M 2 = a ~ asa ~ 4M/j/m + 
C(/i 3 ), corresponding to the superradiance threshold [Eq. <j6j) ] 
when ujr ~ fi. For polar modes we show two different fitting 
functions. Fit I corresponds to Eq. (O with 7_n ~ 20 ± 10, 
i.e, Mlji = 20(a - 2r+u){M^) 7 . Fit II is given by Mcui ~ 
(a - a SR ) MM/x) K ° + r ? ia(M/i) K1 ] with « -6.5 ± 2, 771 « 
2.1 ± I, K « 6.0±0.I, ki « 5.0±0.3. While fit I is physically 
more appealing [26|], fit II does a better job at reproducing 
our numerical data in the whole instability region. 



In Fig. [2] we compare the axial and polar instability 
windows for m v = 10 _20 eV, including our estimated er- 
rors. In the polar case we use two different fitting func- 
tions, to bracket uncertainties. The polar instability win- 
dow clearly depends on the chosen fitting function, but 
some general conclusion can be drawn: (i) axial bounds 
on m v are typically less stringent than polar bounds by 
up to one order of magnitude or more; (ii) different fitting 
functions translate into bounds on m v that differ by a fac- 
tor of a few; (iii) the instability windows in the polar case 
extend down to a ~ 0.1. Thus, essentially any supermas- 
sive BH spin measurement would exclude a considerable 
range of m v ; (iv) from Eq. © it is straighforward to ob- 
tain a conservative upper bound for the boson mass that 
can be excluded by a given observation: 



mah 
4M 



3.34 x 10 



. 19 . 1O*M 
ma — — — ev . 



M 



(8) 



This bound is valid independently of uncertainties in the 
fitting functions. 

Thus, existing measurements of supcrmassivc BH spins 
rule out vector field masses in the whole range 10 -20 cV < 
m„ < 10~ 17 eV. The best bound comes from Fairall 9 
[38| . for which the axial instability implies a conserva- 
tive bound (including measurement and numerical er- 
rors) m v < 4 x 10 -20 eV when we compare the instability 
timescale to the Salpeter time, and m v < 2 x 10 -20 eV 
if we do not consider accretion. This result is of great 
significance, since it is two orders of magnitude more 
stringent than the current best bound on the photon 



mass, m 7 < 10 18 eV [39J. If the largest known su- 
permassive BHs with M ~ 2 x 10 10 M Q [H EJ were 
confirmed to have nonzero spin, we could get bounds as 
low as m v < 10~ 22 eV. 

Nonlinear effects, other couplings. An important 
ingredient that was not taken into account in our study is 
the nonlinear evolution of the instability, that can mod- 
ify the background geometry. Photon self-interactions 
are very weak, being suppressed by the mass of the elec- 
tron. Therefore, it is quite likely that the outcome of 
the instability will be a slow and gradual drainage of 
the hole's rotational energy. Another important issue 
is whether the coupling of accreting matter to massive 
bosons can quench the instability. In principle massive 
photons (unlike hidden U(l) fields, for which the inter- 
action with matter is very small) can couple strongly to 
matter. However it is unlikely that this will significantly 
affect the superradiant instability discussed here, for two 
reasons: (i) the unstable modes are large-scale coherent 
modes whose Compton wavelength is of order the BH size 
or larger, and accretion disks are typically charge-neutral 
over these lengthscales, so any possible coupling with or- 
dinary neutral matter is incoherent and most likely in- 
efficient; (ii) accretion disks are often localized in the 
equatorial plane, and therefore they can affect at most 
some (but not all) unstable modes. The investigation 
of the superradiant instability in the presence of matter 
requires further work, but these arguments suggest that 
estimates in vacuum should be robust. Spin measure- 
ments for slowly accreting BHs (such as the BH at the 
center of our own galaxy) are presumably the most reli- 
able, being rather insensitive to details of the interaction 
of vector fields with matter. 

Conclusions. The results discussed here show that BHs 
offer the exciting possibility to constrain particle physics 
and to set stringent upper bounds on the mass of bosonic 
fields. Our method can be generalized to other fields 
and applied to other backgrounds, such as Kerr-Newman 
BHs or higher-dimensional BHs in TeV-gravity scenar- 
ios. Numerical simulations of either linearized or non- 
linear perturbations of the Kerr geometry are necessary. 
In numerical studies of superradiant instabilities a linear 
perturbation analysis proved to be very useful, e.g. for 
choosing suitable initial data (22[. Numerical methods 
may give us a better understanding of the maximum in- 
stability rate, of the nonlinear evolution and of the end 
state of the instability. These are clearly important top- 
ics: for example, a more accurate quantitative analysis of 
the polar sector could improve the bounds on the photon 
mass by by up to one order of magnitude or more. 
Acknowledgments. We wish to thank the Ax- 
iverse Project members (especially Hideo Kodama) 
and Antonino Flachi for valuable discussions, and 
the referees for their useful suggestions. This work 
was supported by the DyBHo-256667 ERC Start- 
ing Grant, the NRHEP-295189 FP7-PEOPLE-2011- 
IRSES Grant, and by FCT - Portugal through PTDC 
projects FIS/098025/2008, FIS/098032/2008, CTE- 



5 



ST/098034/2008, CERN/FP/123593/2011. A.I. was 
supported by JSPS Grant-in- Aid for Scientific Research 
Fund No. 22540299 and No. 22244030. E.B. was sup- 
ported by NSF Grant No. PHY-0900735 and NSF CA- 
REER Grant No. PHY-I055I03. P.P. acknowledges 
financial support provided by the European Commu- 
nity through the Intra-Europcan Marie Curie contract 



aStronGR-2011-298297 and the kind hospitality of the 
University of Rome "Sapienza" and of the International 
School for Advanced Studies (SISSA) in Trieste. V.C. 
thanks the Yukawa Institute for Theoretical Physics at 
Kyoto University for hospitality during the YITP-T-11- 
08 workshop on "Recent advances in numerical and ana- 
lytical methods for black hole dynamics" . 



and J. March-Russell, |Phys.Rt|v. D81, 123530 (2010) 


arXiv:0905.4720 [hep-th] 






A. Arvanitaki and S. Dubovsky 


Phys.Rev. D83, 044026 


(2011),|arXiv:1004.3558 [hep-th 





939 (2010), arXiv:0809.1003 [hep-ph] 

M. Goodsell, J. Jaeckel, J. Redondo, and A. Ringwald, 

|JHEP 0911, 027 (2009)| |arXiv:b909.0515 [ he P-phH 

J. Jaeckel and A. Ringwald, Ann. Rev. Nucl. Part. Sci. 60, 

405 (2010) 1arXivTT0 02.0329 [hej>ph]| 

P. G. Camara, L. E. Ibanez, ar}d F. Ma rchesano, JHEP 

1109, 110 (2011) 1arXiv:1106.0(D60 [hep-th] | 

W. H. Press and S. A. Teukolsky, Nature 238, 211 (1972) 

T. Damour, N. Deruelle, and R. Ruffini, Lett.Nuovo 

Cim. 15, 257 (1976). I 

T. Zouros and D. Eardley, | Anna[ls Phys. 118, 139 (1979) 
S. L. Detweiler, |Phys. Rev. D22, 2323 (1980)| 
V. Cardoso , O. J. Dias, J. P. ^e mos^ and S. Yoshida, 
|Phys.Rev. D70, 044039 (2004|)[ arXiv:hep-th/0404096 
[hep-th] 

V. Cardoso and S. Yoshida, Ij l HEP 0507, 



|arXiv:hep-th70 502206 [hep-th] 



S. R. Dolan, |Phys.Rev. 
|arXiv:0705.2880 |gr-qc||~ 



(2005) 



)76, 084001 (2007) 



J. Rosa,|jH EP 100 6, 015 (201Q)| arXiv:0912.1780 [hep 
-th]. 

V. Cardoso, S. Chakrabar ti|, P. Pani, E. Berti 
and L. Gualtieri, [Phys.Rev. Lett. 107, 241101 (2011) 
|arXiv:1109.6021 [gr-qcj| ~ 

S. Hawking, Mon.Not.Roy.Astrln.Soc. 152, 75 (1971). 
Y. B. Zel'Dovich and I. D. Novikov, Astron.Zh. 43, 758 
(1966). 

B. J. Carr and S. W. Hawking, Mon.Not.Roy.Astron.Soc. 
168, 399 (1974). 

N. Yunes, P. Pani, and V. Cardoso, Phys.Rev. D85, 

102003 (2012), |arXiv:1112.3351 lir^qcj " 

J. Alsing, E. Berti, C. M. V fill, and H. Zaglauer, 



|Phys.Rev. D85, 064041 (2012)||arXiv:1112.4903 [gr-qc]| 
H. Kod ama and H. Yoshino, [(2011), arXiv: 1108. 1365 
[hep-th]. 

H. Yos hino and H. Kodama, Pr |og.Theor.Phys. 1 
(2012), 



arXiv: 1203.5070 



gr-qc 



G. Mocanu and D. Grumiller, 



[24] 

[25] 

[26] 

[27] 

[28] 
[29] 
[30] 

[31] 



[32] 
[33] 

[34] 

[35] 

[36] 
[37] 



(2012)1 |arXiv:1203.4681 [astro-ph.CO]| 

D. Gjrt'tsov, G. P omerantseva, and G. Chizhov, 

ISov.PHys.J. 27, 697 (1984)] 

C. Heideiro, M. O. Sampaio, and M. Wang, Phys.Rev. 



[39] 
[40] 

[41] 



arXiv:1110.2485 [gr-qc 



D85, 024005 (2012), 

J. G. bosa and S. R. Dolan, P hys.Rev. ^85, 044043 

(2012) narX~iv:1110.4494 [hep-th] | 

R. rjonoplya, |Phys.Rev. D73^ 024009 



arXiv:gr-qc/0509026 [gr-qc] 
Y. Ko 



(2006) 



Ko 



ima, [ Phys.Rev. D46, 4289 (1992)| 
ima, |Astrophys.J. 414, 247 fl 993) 



Progress of Theoretical Physics 90, 977 



Y 

Y. Kcjima, 
(1993) 

P. Pani, V. Cardoso, L. Gualtieri, E. 
and A. Ishibashi, (2012), Codes available 



Berti, 
here: 



http: //blackholes . ist .utl . pt / ?pag e=Files) 
http: //www.phy . olemiss . edu/~berti/qnms ,html| 

arXiv: 1 209.0773 [gr-qc] | " 

S. Chandrasekhar and V. Ferrari, Proc.Roy.Soc.Lond. 
A433, 423 (1991). 

E. B^rti, V. Cardoso, and A. O. Starinets, 
| Class.( juant.Grav. 26, 163001 (2009)[ arXiv:0905.2975 



[gr-qc 

A. A. Starobinsky, Zh. Eksp. Teor. Fiz 64, 48 (1973), sov 
Phys. JETP 37, 28 (1973) 

E. Bern and M. Volonteri,|Astrophys.J. 684, 822 (2008) 



arXiv:0802.00 25 [astro-ph 
J. M. Bardeen, Nature 226, 64 (1970). 
L. Brenneman, C. Reynolds, M. Nowak, 
M. Tippe 



et al. 



103 



R. Reis, 
(2011)1 



|Astrophys.J. 736, 
arXiv:jL104.1172 [astro-ph.HE[| 

K Scjimoll, J. Mill er, M. Volonteri, E. Cackett, 
C. Reynolds, et al, |Astrophys.J. 703, 2171 (2009)| 
arXiv:|)908.0013 [astro-ph.HEj| 

J. Beringer et al. (Particle Data Group), Phys.Rev. D86, 
01000^ (2012) - 

N. J. McConnell, C.-P. Ma, K. Gebhardt, S. A. Wright, 
J. D. Murphy, T. R. Lauer, J. R. Graham, and 
D 



O. 



Richstone, [Nature (London) 480, 215 (2011) 
|arXiv: | L112.1078 [astro-ph. CO j| 
McConnell 



N. J 
hardt, 



Phys.Rev. D85, 105022 



C.-P. Ma, J. D. Murphy, K Geb- 
T. R. Lauer, et a/,.,|Astrophys.J. 756, 179 (2012")] 



arXiv: [20Tl620 [astro-ph.COj 



